Diffusion of Asymmetric Swimmers 
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Particles moving along curved trajectories will diffuse if the curvature fluctuates sufficiently in 
either magnitude or orientation. We consider particles moving at a constant speed with either 
a fixed or with a Gaussian distributed curvature magnitude. At small speeds the diffusivity is 
independent of the speed. At larger particle speeds, the diffusivity depends on the speed through 
a novel exponent. We apply our results to intracellular transport of vesicles. In sharp contrast to 
thermal diffusion, the effective diffusivity increases with vesicle size and so may provide an effective 
means of intracellular transport. 
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The thermal Stokes-Einstein diffusivity of a sphere de- 
creases as the particle radius R increases Q. For this 
reason, while diffusive transport is used for individual 
molecules within living cells Q, larger objects such as 
vesicles and pathogens often use active means of trans- 
port. While many intracellular vesicles appear to be 
transported by molecular motors directed along existing 
cytoskeletal tracks 0, , undirected actin-polymerization 
mediated vesicle transport has been reported in some 
endosomes, lysosomes, other endogenous vesicles, and 
phagosomes [3, ■ Active transport is also observed in 
the actin-polymerization-ratchet motility of certain bac- 
teria [J, |l4j and virus particles |fj within host cells. It 
is important to characterize the transport properties of 
vesicles that are not moving along pre-existing cytoskele- 
tal tracks. 

Existing discussions of the motion of actively propelled 
microscopic particles, or "swimmers" , assumes that in 
the absence of thermal fluctuations particles would move 
in straight trajectories 0, Q- Thermal rotational diffu- 
sion will then randomly re-orient the trajectory [lj, so 
that over long times diffusive transport will be observed. 
However in actin-polymerization based motility, particles 
appear to be attached to their long actin tails |8| which in 
turn are embedded in the cytoskeleton [9|. While ther- 
mal fluctuations will thereby be severely reduced, the 
actin-polymerization itself is a stochastic process with 
its own fluctuations 0, [^} . These intrinsic fluctuations 
can explain the observed curved trajectories, as well as 
the variation of the curvature over time 13] . The dif- 
fusivity of such asymmetrically moving particles has not 
been previously explored. 

In this letter, we study asymmetric swimmers that 
would move in perfect circles in the absence of fluctua- 
tions. We examine both a "broken swimmer" with a fixed 
curvature magnitude and an axis of curvature that is re- 
oriented by fluctuations (rotating curvature, RC), and a 
"microscopic swimmer" with a normally-distributed cur- 
vature that is spontaneously generated by fluctuations 
(Gaussian curvature, GC). In both of these systems, fluc- 



tuations lead to diffusion at long times. We use computer 
simulations to measure the diffusivity of these systems 
as a function of the root-mean-squared curvature Ko, 
the particle speed v, and the timescale characterizing the 
curvature dynamics r. 

We obtain some exact results from polymer systems, 
where each polymer configuration represents a possible 
particle trajectory. Indeed, a broken swimmer with a 
fixed curvature magnitude in d = 3 is exactly analogous 
to the hindered jointed chain discussed by Flory and 
we thereby recover the entire scaling function exactly. In 
that case, the diffusivity is independent of particle speed 
v. For Gaussian curvatures and for systems in restricted 
geometries (d — 2), the polymer analogy gives us the dif- 
fusivity only in the limit of slow speeds. At larger speeds, 
our simulations over 5 decades of speed show that diffu- 
sivity depends on particle speed with a non-trivial expo- 
nent A. The diffusivity appears to be dominated by the 
occasional long straight segments of trajectory that occur 
when the curvature is small. Scaling arguments based on 
this observation are consistent with the measured expo- 
nent \2d = 0.98 ± 0.02 in d = 2, but do not recover our 
measured exponent A3d = 0.71 ± 0.01 in d = 3. 

A curved path has a curvature magnitude K = 1/R, 
where R is the instantaneous radius of curvature. If 
we describe a particle trajectory by a position r(t), 
then the vector curvature is defined by the cross-product 
K = v x v/w 3 , where v = |v| is the speed and the dot ' 
indicates a time-derivative. For uniform motion around 
a circle, R is the radius of the circle, and K is oriented 
perpendicular to the circle along the axis. We consider 
particles moving at a constant speed and with an instan- 
taneous curvature K, so that r = v and v = —w x K. 
For "rotating curvature" dynamics (RC) we fix the cur- 
vature magnitude |K| = Kq but allow the curvature to 
randomly rotate around the direction of motion: 



K 



RC 



£vxK 



(1) 



where the unit-vector v = v/v, the Gaussian noise £ has 
zero mean, and (£(*)£(*')) = 28(t — t'j/r with a char- 
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acteristic timescale r. This represents the simplest de- 
scription of a mesoscopic swimmer that has a "locked-in" 
curvature due to, e.g., an asymmetric shape. For "Gaus- 
sian curvature" dynamics (GC) the curvature magnitude 
changes as well: 

K GC = -K/r + £ (2) 

where the noise £ is perpendicular to v with zero mean 
and (£(t) • £(*')) = S(t - t')K^/r, such that <K 2 ) = K 2 . 
This represents the simplest description of a microscopic 
swimmer "trying" to swim in a straight line subject to 
intrinsic fluctuations in the motion. The resulting curva- 
tures are Gaussian distributed in each component. For 
particles restricted to two-dimensions with either RC or 
GC dynamics, we only use the normal (z) component of 
the vector-curvature to update the velocity within the 
plane, i.e. v = — vv x zK z in d = 2. 

There are two natural timescales. We explicitly in- 
troduce r, which controls the noise correlation and so 
sets the timescale over which the curvature changes. 
There is also the inverse of the angular rotation rate, 
t c = l/(vK ). Diffusion will only be observed for elapsed 
times t much greater than any other timescale in the sys- 
tem, i.e. t 3> t c and f > t. The diffusivity of a particle 
is given by D = (r 2 ^ /(2dt) in the limit as the elapsed 
time t — > oo, in spatial dimension d. 

A polymer chain with fixed bond lengths (£) and an- 
gles (Of), and with independent bond rotation potentials 
(V(4>f)) UIl is statistically identical to the continuous 
RC trajectory in 3d if for a discrete time-step At we take 
i = vAt. The end-to-end distance for a long n-bond 
polymer is (V 2 ) = n£ 2 C n . The correspondence is com- 
plete as the elapsed time t = nAt — > oo. The bond 
and dihedral angles determine = (1 + cos#/)(l + 
(cos0/))/[(l — cos#/)(l — (cos0/))] 1 lj . Swimmers fol- 
low continuous paths, so we take the limit of small At 
and fix the polymer rotation angle from the curvature in 
that limit by 9f — K^vAt, and rotate the curvature by 
< <fi 2 >= 2 At It in agreement with Eqn. ^ In the limit 
At — > we recover the exact result D = 1/(3Kqt) in 
d = 3. Remarkably, D is independent of v. 

For a GC trajectory in d = 3, there is no obvious 
polymer analogy since the curvature magnitude evolves 
with time. In the limit of r — > however, the curvature is 
independently Gaussian distributed at every point along 
the trajectory and the diffusivity can be extracted from 
the "worm-like chain" polymer model originally solved 
by Kratky and Porod O- The result is D = 1/(3J^t) 
in the limit of small r. Note that the diffusivity diverges 
as 1/r so this is the leading asymptotic dependence for 
small t. We obtain the same diffusivity for both RC and 
GC dynamics for small r. 

We use these exact results to define natural dimension- 
less scaling functions for the diffusivity of microscopic 



swimmers: 

D R d,Gd{v) = DK 2 t (3) 

where the index Rd or Gd indicates both the dynamics 
(RC or GC) and the spatial dimensionality d, and 

v = vKqt (4) 

is a dimensionless speed. In terms of these scaling 
functions we have Drz{v) = Dg3(0) = 1/3. The 
same Kratky-Porod approach in d = 2 gives Dr2 (0) = 
D G 2(0) = 1. 




FIG. 1: a) Particle trajectory with GC dynamics in d = 2 
with v — 0.1. The particle does not complete a circular loop 
before K changes significantly, b) With v = 100. The particle 
can complete many circular loop before K changes, however 
straight segments are seen when |K| is small. The result is a 
characteristic "knotty wool" appearance. In both cases Ko = 
1. 

We have simulated the trajectories of large-numbers of 
independent particles with RC and with GC dynamics. 
For fixed v and Kq, we varied r to explore the scaled 
velocity v = vKqt over 5 orders of magnitude. For each 
v, we averaged over the trajectories of at least 1000 par- 
ticles. We explicitly integrated the dynamical equations 
using a simple Euler update with a small timestep At. In 
all cases t ^> t ^> At and t S> t c S> At, with separation 
of timescales by factors of 10 — 100. Systematic errors 
due to At and t are below our noise levels, and statistical 
errors (when not shown) are smaller than the size of our 
plotted points. Consistently, our numerical results agree 
with all exact results from the polymer analogy. We illus- 
trate the trajectories that we observe in d = 2 in Fig. ^ 
with both small and large scaled speeds v. In both cases 
the curvature Kq = 1, but particles only complete loops 
at large v. Qualitatively similar trajectories are seen in 
d = 3 with GC curvature dynamics. 

In d = 2, shown in Fig. |2 both rotating curvature 
(open circles) and Gaussian curvature (filled circles) ap- 
proach their asymptotic value of Dr2{Q) — Dq2{Q) = 1 
at small v. At v w 1 there is a sharp cross-over to a 
large-w power-law regime, characterized by an exponent 
\2d where Dr2 ~ Dgi ~ v X2d for large v. We show the ef- 
fective exponents X e ff = A log (DK 2 t)J A log v between 
consecutive points in the inset of Fig. |2 as well as the 
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FIG. 2: Dimensionless diffusivities D = DKqt for rotating 
(open circles, Dr2) and Gaussian (filled circles, Dg2) curva- 
ture dynamics in d — 2, plotted against dimensionless par- 
ticle speed v = vKor. Also shown with solid lines are the 
small v asymptote D = 1 and the large v best-fit asymptote 
D ~ v X2d . The inset shows the effective exponents, with the 
solid line indicating the best-fit \2d = 0.98 ± 0.02. 
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FIG. 3: Dimensionless diffusivities D = DKqt for rotating 
(open circles, Dm) and Gaussian (filled circles, Dg3) curva- 
ture dynamics in d = 3, plotted against dimensionless particle 
speed v = vKqt. Solid lines show the exact result Drz = 1/3, 
as well as the large v power-law asymptote Dei ~ v Xad . 
The inset shows effective exponents between sequential points, 
with a solid line indicating the best-fit \$d = 0.71 ± 0.01. 



best fit exponent \2d = 0.98 ± 0.02. We fit \2d from the 
large-i; GC data only, due to the systematic cross-over 
remaining in the RC data even at large v. 

Simulations in d — 3 with rotating curvature (RC) dy- 
namics leads to a diffusivity in excellent agreement with 
the exact result from polymer physics, = 1/3, as 

shown by open circles in Fig. [3] Gaussian curvature dy- 
namics (Dg3, filled circles) has the same behavior for 
small v, but exhibits a sharp crossover at 5 w 1 to a 
power-law regime Dqz ~ v X3d for large v. We find the 
best-fit exponent is \$d = 0.71 ± 0.01, as shown by the 
solid line in the inset of Fig- EI Because \3d < 1, this scal- 
ing curve may be used to uniquely identify the dynamical 
timescale r if D, Kq, and v are measured experimentally. 

Is there a simple way of understanding the asymptotic 
behavior of D1 For RC dynamics in d = 3 the instanta- 
neous curvature does not change in magnitude even while 
the curvature axis wanders. The particle will go in a cir- 
cular trajectory, not contributing to diffusivity, until the 
curvature axis wanders significantly. The result is a ran- 
dom walk with step size given by the radius of curvature 
Ar ~ l/-Ko and an interval between steps of r, leading to 
D ~ 1/(Kqt). This qualitatively explains why the exact 
result Drt, = 1/3 is independent of v. 

It is more difficult to understand the D ~ v x behavior 
for large v in the other systems. We start with a simple 
scaling argument based on the assumption that the rela- 
tively straight segments shown in Fig.^b) dominate the 
diffusivity. The interval between periods of small curva- 
ture should be on order the autocorrelation time r. The 
length Ar of the straight segments are determined by 
how long the interval of small curvature lasts, At, since 



Ar w vAt. For the segment to be straight, the curvature 
must be less than the inverse length, i.e. K max < 1/Ar. 
The fraction of the time we have small curvature below 
K m ax in magnitude should be proportional to the prob- 
ability of having curvature below K max . In d = 2 only 
the normal component of curvature affects the dynamics, 
so that P(K) « const for K <C K$. This applies both 
to GC and RC. We therefore expect At ~ TK max / Kq. 
We maximize K max to maximize the contribution to 
D ~ Ar 2 /r and find Dq2 ~ Dr2 ~ v as v — > 00. 
This indicates that A2<j = 1, which is consistent with 
our best-fit value \2d = 0.98 ± 0.02. However, in d = 3 
for GC dynamics the same argument leads to X^d = 2/3 
since two Gaussian distributed components of the cur- 
vature gives P<(K m?x ) = J*"" dKP(K) ~ K max /K$ 
for K max <gC Kq. This is inconsistent with our measured 
value of \3d = 0.71±0.01, with a significant 4er variation. 

At what radius R c does a small spherical particle 
achieve a higher diffusivity by actively swimming, as 
compared to passive thermal diffusion characterized by 
D T = k B T/(6nriR) Q? We can answer this question 
within the context of actin-polymerization based motil- 
ity of small intracellular particles, since the size depen- 
dence of Kq, v, and r is known, at least approximately. 
With the approximation that n propulsive actin filaments 
are randomly distributed over a particle of size R, the 
curvature of the trajectory will be Kq oc l/(Ry/n) [l3| . 
With a size-independent surface-density of filaments we 
obtain Kq w A/R 2 , with a constant of proportionality A. 
By observations of Listeria monocytogenes we estimate 
A w 0.1[im 0. We also conservatively assume size- 
independent values for cytoplasmic viscosity n ps 3Pa • s 
[H fl3l|. speed v ps 0.1/im/s, and autocorrelation decay 
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time t w 100s |15(. We find that the micron-scale bac- 
terium L. monocytogenes has v ~ 1, so that smaller par- 
ticles will have v > 1. Using the large z) asymptotic 
behavior of D G3 shown in Fig. El D w 0.41u A3d /(^r), 
and the size-dependence l<"o ~ A/R 2 , we obtain 

L» G3 - i?4-2A 3du A 3<J /(^-Aa^l-Asd^ ( 5 ) 

with a measured = 0.71 ± 0.01. In dramatic contrast 
to thermal diffusion, D increases with increasing parti- 
cle size. Comparing with Dt we find that for all sizes 
above R c « 80nm a particle will have a higher diffusivity 
by actively swimming by the actin-polymerization mech- 
anism than by passive thermal diffusion. Provocatively, 
this is in the middle of the vesicle-size distribution seen 
in neural systems [lfil ]. 

Our treatment of microscopic swimmers has ignored 
thermal fluctuations. A "rocket" traveling straight at 
speed v that is re-oriented only by thermal effects will 
have D u — 4irr]R 3 v 2 / (SksT) [lj. In comparison with 
our results for D, we find that D < D u for particles 
larger than R u « 0.07nm. For actin-polymerization 
based motility, intrinsic fluctuations appear to dominate 
thermal fluctuations at the particle sizes where active 
transport is advantageous. 

In summary, we find that diffusivities of asymmetric 
microscopic swimmers depends on whether the swimmers 
are restricted to 2d or 3d, and whether they have fixed 
asymmetries (RC) or the asymmetries are spontaneously 
generated (GC). Diffusivities are independent of particle 
speed at low speeds, in agreement with analogous poly- 
mer systems. At higher speeds an anomalously large dif- 
fusivity is observed that depends on the particle speed 
by v x where \ 2 d = 0.98 ± 0.02, in agreement with a scal- 
ing argument for X 2 d = 1. However X S d = 0.71 ± 0.01, 
which significantly differs from our simple scaling result 
in d = 3. We apply our results to intracellular bac- 
teria, virus particles, and vesicles that move via actin- 
polymerization. We find that diffusivities due to asym- 
metric swimming exceed thermal diffusivities for particles 
larger than approximately 80nm. As a result asymmet- 
ric swimming may provide a viable intracellular trans- 
port mechanism even for vesicle-sized particles. We find 



that for the relevant dynamics (GC in d = 3), diffusivi- 
ties should increase with particle size, speed, and filament 
turnover rate, and also with smaller curvatures for a given 
size. It is interesting that the bacterium Rickettsiae rick- 
ettsii exhibits actin-polymerization intracellular motility 
with smaller intra-cellular speeds but straighter trajecto- 
ries 0, ^3 — raising the question of whether maximal 
diffusivity is selected for in this or other biological sys- 
tems. 
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